Role of c-axis pairs in V2O3 from the band-structure point of view. 
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The common interpretation of the LDA band structure of V2O3 is that the apparent splitting of 
the ai 3 band into a low intensity structure deep below the Fermi energy and a high intensity feature 
above it, is due to the bonding-antibonding coupling of the vertical V-V pair. Using tight-binding 
fitting to -as well as first-principles NMTO downfolding of- the spin-up LDA+U ai g band, we show 
that there are other hopping integrals which are equally important for the band shape as the integral 
for hopping between the partners of the pair. 



A few years ago, Park et al. reexamined the socalled 
spin 1/2 model for V2O3 yj. Based on polarization 
dependent x-ray absorption measurements they showed 
that, for all phases, the vanadium 3+ ion (d 2 ) is in the 
spin 1 state. They also demonstrated that this state 
is a mixture of e g e g and e£a lff configurations, with the 
former having the larger weight, especially at low tem- 
peratures. Recall, that the t2 g orbitals, which are pdir 
anti-bonding with the O p orbitals on the surrounding 
octahedron, lie below the pda anti-bonding e g orbitals 
and are split by a trigonal distortion into low-lying, dou- 
bly degenerate orbitals and a higher-lying a\ g orbital. 
The picture presented by Park et al. is hardly consis- 
tent with the classical vertical-pair assumption that the 
bonding-antibonding splitting of the et lff orbitals of the 
V-V pair places the energy of the bonding orbital well 
below that of the orbitals 0, 0] . This is the assump- 
tion which 25 years ago led Castellani et al. to suggest 
the spin 1/2 model where for the V-V pair two electrons 
fill the bonding a\ g orbital and the two remaining elec- 
trons 1 / 4- fill the four orbitals [4j . It is now generally 
recognized that the spin 1/2 model is incorrect. Never- 
theless, the vertical pair remains a popular starting point 
for current attempts to calculate the electronic structure 
of V2O3 H,IEI3- A comprehensive review of the latest 
experimental and theoretical results in this field can be 
found in a recent paper by Di Matteo, Perkins and Natoli 

In the present paper we study the dispersion of the 
a\ g band obtained from a modern LDA+U calculation by 
performing a tight-binding analysis. Our motivation for 
doing this is to obtain information regarding the small- 
est cluster that one can use in the model calculations 
while still preserving the most important aspects of the 
band structure. In spite of the fact that some literature 
exists providing the qualitative hint concerning this, we 
considered it to be rather important to check this with 
more recent methodology (e.g. the downfolding tech- 



nique) providing the quantitative estimates. Note that, 
although the spatial orientation of an a\ g orbital is ac- 
tually in favour of a molecular-like picture, our analysis 
shows that there are other hopping integrals which are 
equally important for description of the main features of 
the band structure as the hopping integral between the 
partners of the c axis pair (ti). Though these inter-pair 
hopping integrals are smaller than ti, the contribution is 
proportional to the number of neighbours and that makes 
them rather significant. 

Since we are interested in understanding the relative 
importance of the hopping between pairs and non-pairs, 
we consider the band structure of ferromagnetic V2O3 
in the high-temperature corundum structure. Although 
this phase does not exist in nature, it can provide a good 
estimate of an upper bound for the hopping integrals in 
this compound for the following reasons: First of all, 
comparing with antiferromagnetic phases, the ferromag- 
netic one has the largest band-width Secondly, the 
distance between the partners of a vertical pair is shorter 
in the corundum structure than in the low-temperature 
monoclinic phase [Tl| . Therefore the hopping integral 
between partners of the pair in this structure should be 
maximal. 

In FigQ]we show the LDA and LDA+U spin-up band 
structure of V2O3 in the energy range of the 12 ti g bands 
(4 V atoms per cell). Comparing these two pictures one 
can see that the LDA+U band structure calculated for 
U=3eV and J=0.8eV is essentially a rigid shift of the 
band down in energy and a\ g up so that the former is 
almost completely below the Fermi energy and the later 
is above it. However, we note that as a result the mixing 
between e g and a\ g bands is suppressed. Nevertheless, 
it is clear from FigQ] that this mixing does not come 
from the hybridization between different orbitals of the 
atoms in the c-axis pair. This warns us already about 
the importance of the neighbours other than the partner 
in the pair. 
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FIG. 1: LDA (upper panel) and LDA+U (lower panel) spin- 
up band structure of ferromagnetic V2O3 in the corundum 
structure with 2 formula units per cell. The amount of a\ g 
character is indicated by the width of so-called fat bands (dot- 
lines). We note that T-Z is along the direction of the verti- 
cal pair (see Fig. 3.11 in 9] and Fig.l in [l|]). r=(0,0,0); 
£=(1/2,1/2,1/2); L=(0,l/2,0); F=(l/2,l/2,0). The zero of 
energy is at Fermi energy. 



The LDA+U band structure, which yields a spin 1 
state, has the advantage that, as already mentioned, 
the empty a\ g band is practically separated from the 
full band. Of course, this depends on the values of 
the parameters U and J used in the calculations. Ac- 
cording to Solovyev et al. the calculated value of the 
screened parameter U and Hund's rule exchange J for 
V 3+ ion in LaVC>3 are 3eV and 0.93eV, respectively 
[l2|. On the other hand, an empirical estimate by Marel 
and Sawatzky, based on gas-phase multiplet splittings of 
the 3d series, shows that in the case of V 3d J is about 
O.74eV0^ In the present work we use J=0.8eV as es- 
timated by Tanabe and Sugano for the free V 3+ ion 
00 

Before we discuss our results we note that, the width 
of the aig band reaches its maximum at T, where it is 
about 2eV. One might believe that this width is mainly 
caused by the bonding-antibonding interaction between 
the vanadium pairs along the c-axis, and that the inter- 
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FIG. 2: Vanadium neighbors in corundum structure which are 
taken into account in the tight-binding model (a). Definition 
of hopping parameters (b) ii between the pair VI and V3, 
(c) fa between VI and V2, (d) £3 between VI and V3' and 
<4 between VI and V4. The small spheres denote the oxygen 
ions to illustrate the importance of the bond angles. 



action between pairs is weak. In the simplest such pic- 
ture only two hopping parameters would be important: 
A large intra-pair hopping parameter, which should give 
the most of the bandwidth, and a smaller inter-pair hop- 
ping. This would result in a situation where the splitting 
of the aig band at the T-point is primarily determined 
by the value of intra-pair hopping integral. However, one 
notices that the band splits into four almost equally sep- 
arated levels. 

To shed more light on this issue, we carry out a tight- 
binding model calculation where the hopping integrals to 
the first (ti), second (£2), third (£3) and fourth (£4) near- 
est V neighbors are taken into account (Fig|2J. Again, 
t\ is the hopping integral between the atoms of the c- 
axis pair. In Fig. 3 (a) to (d) we demonstrate how each of 
these parameters contribute in the dispersion of a\ g band. 
Namely, switching on merely the hopping parameter t\ 
splits the atomic a\ g level into two doubly degenerate flat 
bands with energy difference 2t\. There is no dispersion 
because the atoms in-between the pairs in the structure 
are missing. In fact, only because these atoms are miss- 
ing, does one see pairs at all. Now the main dispersion 
is caused by the hopping parameter i 3 which yields a 
maximal splitting of at the T-point (the number of 
neighbors which an electron can hop to with £3 is equal to 
3). The hopping parameter t± lifts the degeneracy of each 
of these doubly degenerate bands. Inclusion of t% makes 
the band asymmetric with respect to the position of the 
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FIG. 3: Tight-binding ai g bands. The zero of energy is at the center of the gap between the split bands, (a) ti=-0.25eV, 
t2=0, £3=0, £4=0. The splitting is 2 x ti and there is no dispersion, (b) ii=-0.25eV, i2=0, i3=-0.15eV, t4=0. Note that the 
large splitting at F is dominated by ts and not t-z. (c) ti=-0.25eV, £2=0, t3=-0.15eV, t4=-0.06eV. A small £4, but with 6 
nearest neighbors, is sufficient to cause the large splitting into 4 bands seen in the LDA+U band structure, (d) ii=-0.25eV, 
£2=-0.03eV, t3=-0.15eV, i4=-0.06eV. To obtain the very evident asymmetry between the upper and lower pairs of bands we 
need a small £2- This is now close to the LDA+U band structure. 



initial a\ g level. The final result in Fig. 3(d) looks very 
much like the LDA+U a\ g band. The parameters used 
are: ii=-0.25eV, t 2 =-0.03eV, i 3 =-0.15eV, t 4 =-0.06eV. 
Although t\ is indeed larger than other hopping integrals, 
its influence on the bandwidth is not that large, because 
there is only one nearest neighbour, as compared to 3 for 
t 2 , 3 for £ 3 and 6 for t\. 

It is important to note that these values for the hop- 
ping parameters are not unique. For example, t\=- 
0.5eV, i 2 =-0.03eV, t 3 =-0.1eV and i 4 =-0.04eV would 
also give a small splitting at the Z-point, and large one 
at r or L, as is shown in Fig. 4. This has to do with the 
symmetry of these points. At the Z-point, for instance, 
the splitting between the upper and lower components of 
the a\ g band is 2|3*3 — ty\, for £2=0. Hence, for any given 
value of £3 there are always two values of t\ ~ 3^3 ± <5z/2 
which give exactly the same splitting. Note that £4 does 
not influence the energies at Z. At the L-point, on the 
other hand, the splitting is determined primarily by the 
sum of t\ and £3. Therefore, no set of parameters with 



t\ greater than 0.5eV can reproduce the LDA+U band 
structure at this point (unless £3 and t\ have the opposite 
signs, which is " unphysical" ) . Comparing figures 1 and 4, 
one can easily see that the two sets of parameters, which 
give the same splitting at Z, will give different levels at 
F, and that only the set with t\ =-0.25 eV reproduces 
the accidental degeneray of the two middle levels in the 
LDA+U. However, to reproduce the LDA+U band in 
such detail may not be meaningful as long as all hopping 
integrals beyond £4 are neglected. 

Although our tight-binding study has clearly demon- 
strated the role of the various hopping integrals for the 
dispersion of the a\ g band, and although we can find pa- 
rameters which fit the LDA+U a\ g band, it is difficult to 
select one set of hopping integrals because, to the accu- 
racy expected for our model, different sets can do this. 
The most straightforward way to resolve this problem is 
to use the downfolding procedure of Andersen et. al. 0] . 
This procedure relies on keeping in the first-principles 
NMTO band-structure calculation only the relevant de- 
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FIG. 4: Tight-binding oi 9 bands for ti=-0.5eV, <2=-0 
t 3 =-0.1eV and i 4 =-0.04eV. 
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grees of freedom, in this case the a\ g Wannier-like orbitals 
whose LDA+U spin-up bands lie in the energy range from 
the Fermi level to nearly 2 eV above, and integrating out 
the other degrees of freedom. This naturally takes into 
account re-normalization effects due to the integrated-out 
orbitals. Fourier transform of this few-orbital downfolded 
and symmetrically orthonormalized NMTO Hamiltonian 
provides the hopping matrix elements of the correspond- 
ing tight-binding Hamiltonian. This method provides a 
way of generating Wannier-like functions and their single- 
particle Hamiltonian without any fitting procedure. The 
detailed discussion of such calculations for V2O3, as well 
as comparisons with Hamiltonians proposed previously, 
will be presented elsewhere |l7| . Here we only mention a 
technical point specific to the present application: Since 
the downfolding procedure takes place at a more basic 
level than where U is "added" to the LDA, we need to 
construct the potential which yields the spin-up LDA+U 
band structure. That potential we obtained from the 
LDA potential by shifting its logarithmic-derivative func- 
tions, tp' Rlm (e, s) /ipmm (e> s) , m energy so as to repro- 
duce the spin-up LDA+U band structure. 

From this NMTO downfolding calculation we ob- 
tained the following hopping integrals: ti=-0.30eV, ^2=~ 
0.02eV, i 3 =-0.11eV and £ 4 =-0.05eV, which are close, 
although not identical to those used in Fig. 3(d). As 
Fig. 5(a) shows, the band structure obtained from these 
hopping integrals differs a bit from that in Fig. 3(d), and 
from the upper four LDA+U bands in Fig^ The reason 
is simply that the downfolded band structure shown in 
Fig. 5(b), obtained by downfolding to the the four a\ g 
bands, cannot be reproduced completely with merely 
tij ^2, ^3, and t±. Its Hamiltonian has also non-zero 
higher Fourier components, which is hardly surprising. 
The downfolding calculation thus confirms the gross val- 
ues of the hopping integrals found by tight-binding fitting 
to the first 4 shells, but also points to the need for in- 
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FIG. 5: Tight-binding a\ g bands calculated with ti=-0.3eV, 
t 2 =-0.02eV, t 3 =-0.11eV and £ 4 =-0.05eV (upper panel) and 
those obtained from the downfolding procedure (lower panel). 
Note that the zero of energy in the upper panel corresponds 
to the energy of atomic a% g level but it is at Fermi energy in 
the lower panel. 



eluding longer ranged hoppings to reproduce the details. 

One should keep in mind that the hopping integrals 
discussed above are for vanadium-centered Wannier-like 
orbitals and thus different from the Slater- Koster hop- 
ping integrals for atomic oxygen and vanadium orbitals 
obtained by Mattheiss [l^- Following Harrison one can 
show that the hopping integral between ai n atomic or- 
bitals on the vanadium pair is about 0.8eV |ljj, whereas 
the one in which the oxygen degrees of freedom are inte- 
grated out is much less than that. This reduction is due 
to the anti-bonding character of the pdir interaction. 

We have thus demonstrated that, although the inte- 
gral for hopping between the vertical pair is the largest 
hopping integral, it is not the single most important one 
for the aig bandwidth. This is so because the actual hop- 
ping integrals are not only determined by the direct V-V 
hoping but also evolve via intermediate O 2p orbitals. 
The simple picture where only the hopping parameter 
within the c-axis pair is important is not sufficient to de- 
scribe the aig band in V2O3. Our calculations show that 



the hopping integrals between second, third, and fourth 
nearest vanadium neighbors are equally important. 
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